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INFLUENCE  OF  MAGNETIC  SHEAR  ON  THE  LOWER-HYBRID-DRIFT  INSTABILITY 

IN  FINITE  0  PLASMAS 

I.  INTRODUCTION 

The  lower -hybrid -dr if t  instability  is  believed  to  be  an  important 
microinstability,  in  both  space  and  laboratory  plasmas,  because  of  the 
anomalous  transport  properties  associated  with  it.*  This  instability  is 
driven  by  the  diamagnetic  current  in  an  inhomogeneous  plasma  and  is 
attractive  because  (1)  it  can  be  excited  by  modest  density  gradient  (i.e., 

Hi  ^  ^rai^me  t/4  'u  where  Ln  is  the  scale  length  of  the  density  gradient 
and  r  .  is  the  mean  ion  Larraor  radius)  and  (2)  it  is  relatively  insensitive 

LI 

to  the  temperature  ratio  Tg/Tj  (unlike,  say,  the  ion  acoustic 

instability).  The  mode  originally  gained  interest  as  a  means  to  provide 

anomalous  diffusion  of  particles  in  9  pinch  experiments  during  late-time 

sheath  broadening.  *’-*•  ^  Subsequently,  it  has  been  applied  to  other 

laboratory  devices,  such  as  toroidal  reversed  field  pinches’’*6  and  compact 
7  R 

torii,  and  to  magnetospheric  plasmas.  Experimentally,  the  lower-hybrid- 
drift  instability  has  been  directly  observed  in  a  laboratory  plasma^  and 
has  been  suggested  a  mechanism  to  explain  satellite  observations  of 
fluctuating  electric  fields  in  the  earth's  magnetopause*®  and 
magnetotai l. ** 

At  present,  there  is  considerable  interest  in  the  lower-hybrid-drift 

instability  in  regard  to  its  influence  on  the  dynamics  of  reversed  field 

plasmas.  It  is  being  applied  to  the  field  reversed  experiments  at  Los 
U  12 

Alamos  *  and  to  reconnection  processes  In  the  earth's 

magnetosphere.  *  One  of  the  Important  features  of  a  field  reversed 

plasma  (specifically,  one  containing  a  field  null)  is  that  8  varies  over  an 

enormous  range  in  the  reversal  region.  To  accurately  describe  lower- 

hybrid-drift  waves  in  such  a  plasma,  the  analysis  must  include 

2 

finite  8  effects.  The  two  important  finite  8  effects  that  enter  the 
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problem  are  (1)  the  coupling  of  electrostatic  and  electromagnetic 
oscillations  and  (2)  the  orbit  modification  of  the  electrons  due  to  VB. 
Aside  from  finite  8,  another  potentially  important  feature  in  reversed 
field  plasmas  is  magnetic  shear;  the  reversing  field  could  be  sheared  by  a 
component  of  B  which  is  parallel  to  the  plasma  current  that  generates  the 
field  reversal. 

The  initial  study  of  the  influence  of  magnetic  shear  on  the  lower- 

hybrid-drift  instability  considered  electrostatic  waves  in  a  low  8  plasma 

using  analytical  analysis.^  It  was  found  that  magnetic  shear  stabilized 

the  instability  for  L  <  L  (r, . /L  +  L  /r. . )  where  L,  is  the  scale  length 
J  s  n  Li  n  n  Li  s  6 

associated  with  the  magnetic  shear.  A  numerical  analysis  of  this  problem 
was  presented  in  Gladd  et  al.  (  1977). Davidson  et  al.  (1978)^  extended 
previous  analyses  to  finite  8  plasma.  However,  they  considered  the 

limit  *  0  so  that  electron-wave  resonances  could  be  ignored.  The 
important  finite  8  effect  retained  was  the  coupling  of  electrostatic  and 
electromagnetic  oscillations.  Recently,  iluba  et  al.  (1982)^  have 
developed  a  theory  of  the  instability  in  finite  8  and  Tg  plasmas  containing 
magnetic  shear.  However,  the  results  are  restricted  to  the  parameter 
regime  T^  «  T^  and  8g  <<  1.  Thus,  to  date,  a  general  self-consistent 
theory  of  the  lower-hybrid-drift  instability  in  a  finite  8  and  Tg  plasma 
containing  a  sheared  magnetic  field  has  not  been  developed.  The  purpose  oi 
this  paper  is  to  present  such  a  theory. 

The  major  results  of  this  work  are  the  following.  The  finite  8  effect 
arising  from  the  coupling  of  electrostatic  and  electromagnetic 
perturbations  can  have  a  destabi 1 izing  influence  on  the  instability  in  a 
sheared  magnetic  field.  That  is,  as  8  is  increased  in  certian  parameter 
regimes,  the  growth  rate  of  the  instability  increases.  Physically,  this 
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occurs  because  the  fluctuating  electric  field  associated  with  6A„  inhibits 
electron  flow  along  the  magnetic  field  which  reduces  the  rate  at  which 
energy  can  be  convected  away  from  the  localization  region.  The  finite  3 
effect  of  electron  orbit  modification  has  a  stabilizing  influence  on  the 
instability  in  a  sheared  field.  This  is  due  to  a  VB  electron  drift-wave 
resonance  which  is  a  dissipitat ive  effect.  The  key  parameter  which 
dictates  which  finite  3  effect  is  more  important  is  Tg/Tt.  In  the  limit  Te 
«  Tj  the  electromagnetic  coupling  is  dominant,  while  for  T^  the  VB 
electron-wave  resonance  is  more  important. 

The  structure  of  the  paper  is  as  follows.  In  Section  II  we  present 
the  assumptions,  equilibrium,  and  derivation  of  the  equations  which 
describe  the  lower-hybrid-drift  instability  in  a  finite  3  plasma  containing 
magnetic  shear.  In  Section  III,  the  results  are  presented,  both  analytical 
and  numerical.  And  finally,  the  last  section  summarizes  the  results  and 
discusses  the  application  of  this  work  to  reversed  field  plasmas.  Details 
of  the  calculation  are  presented  in  Appendices  A  and  B. 


II.  THEORY 


A.  Assumptions  and  Equilibrium 

We  consider  a  slab  geometry  plasma  which  contains  inhoraogeneties  in 

the  density  and  the  magnetic  field  in  the  x-direction.  The  temperature  is 

assumed  constant  for  simplicity.  Equilibrium  force  balance  on  an  ion  fluid 

element  in  the  x-direction  requires  that  where 

V,.  =  (v?/2rt.)  3  in  n/dx  is  the  ion  diamagnetic  drift  velocity, 
dill 

1/2 

=  (21^/iik)  is  the  ion  thermal  velocity  and  =  e  B^^/m^c  is  the  ion 

Larmor  frequency.  The  ion-diamagnetic  velocity  can  be  related  to  the  mean 

ion  Larmor  radius  and  the  scale  length  of  the  density  gradient  by 

V.^/v^  =  rLi/2Ln  where  rL1  =  v^ft^  and  Ln  =  (3in  n/3x)  1  is  the  density 

gradient  scale  length.  We  consider  magnetized  electrons  while  the  ions  are 

kept  unmagnetized.  This  is  reasonable  in  treating  the  lower-hybrid-drift 

instability  since  we  are  considering  waves  such  that  Slj  «  to  <<  and 

k2r  2  »  1.  We  assume  that  the  plasma  is  weakly  inhomogeneous,  i.e, 

Li 

2  2 

r2  (3£n  n/3x)  «  1  and  r2  (3fcn  B/3x)  «  1.  The  plasma  g  is 

Le  Le 

arbitrary,  l^e  inhomogeneous  ambient  magnetic  field  is  given  by 

B(x)  =  Bq { [ 1  +  (x-XQ)/LB]ez  +  (x-xo)/Lg  ey}  (l) 

in  the  vicinity  of  xQ  (i.e.  (x-xQ)/Ls  «  l  and  (x-x0)/Lg  «  1  where 

L„  «*  OS/Sx)"1,  $  =  tan"1  (B  /B  )  and  LD  =  (3  In  B  /3x_1)).  Thus,  Lc  is 
s  y  z  o  z  ® 

the  scale  length  characterizing  the  magnetic  shear  and  Lg  characterizes  the 

magnetic  field  gradient  scale  length. 

In  the  absence  of  any  field  inhoraogeneities ,  the  field  configuration 

is  B  *  B  e  •  The  plasma  described  above  is  unstable  to  the  kinetic  lower- 
—  o  z 

hybrid-drift  instability  when  1  >  V,./v.  >  (m  /m. )  '  .  0  The  instability 

d  i  1  e  i 
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is  driven  by  by  cross-field  current  and  is  excited  via  an  ion-wave 

resonance  (i.e.,  invers  Landau  damping).  The  waves  are  characterized  at 

maximum  growtli  by  u>r  ~  Vdi  <  “jlh*  Y  *  “r»  kypes“  and  k  •  B  =  0  where 
1  /2 

p  =  (2T./m  )  /SI  .  For  modes  such  that  k  •  B  *  0 

es  i  e  e  — 

(i.e.,  ka  *  0)  electron  Landau  damping  reduces  the  growth  rates  or 

stabilizes  them,  depending  on  the  magnitude  of  k^. 

The  above  description  of  the  plasma  is  significantly  modified  by 

introducing  a  shear  in  the  magnetic  field  as  given  in  Eq.  1.  The  magnetic 

field  rotates  in  the  y-z  plane  (see  Fig.  1)  as  a  function  of  x.  At  xQ  we 

see  that  k^  =  0  while  at  x  »  xj,  k^  *  0  .  Thus  the  dispersive  properties 

of  the  plasma  are  also  a  function  of  x.  We  introduce  the  effect  of  the 

magnetic  shear  (i)  locally  through  k  +  k  (x)  =  k  +  k  (x-x  )/L  and  (ii) 

Ct  z  zu  y  os 

globally  by  replacing  ikx  by  3/3x.  We  note  that  the  magnetic  shear 

distorts  the  particle  orbits  in  a  uniform  magnetic  field  and  introduces  a 
kinematic  drift  term.17  The  orbital  effects  of  shear  will  be  considered 

elsewhere. 
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Figure  1 

Schematic  ot  a  sheared  magnetic  fi 
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B.  Dispersion  Equation 

An  outline  of  the  derivation  of  the  dispersion  equation  which 
describes  the  lower-hybrid-drift  instability  in  a  finite  0  plasma 
containing  a  sheared  magnetic  field  is  presented.  Details  of  the  analysis 
are  presented  in  Appendices  A  and  B.  The  Maxwell  equations  for  the 
perturbed  fields  are 


7  •  6E  *  4ir  6p 


(2) 


7  x  6B 


/  i  asE 

is  «  j  +  ± 

c  -  c  3t 


(3) 


where  6p  and  6J  are  the  perturbed  charge  density  and  current, 
respectively.  Equations  (2)  and  (3)  can  be  rewritten  as 


7  x  7  x  6A  + 


325A 


3t  ‘ 


+  i  7  Mi  =  4i 

c  3t  c 


(4) 

(5) 


where 


6B  =  7  x  6A 

36A 

6E  =  -  76*  -  -  — 

-  y  c  3t 

relate  the  perturbed  fields  to  the  perturbed  potentials,  and 

6n  *  /  d^v  Sf 
o'  o 


(6) 

(7) 


(8) 
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6J  =  e  fd^v  v  5f 
~  a  ‘  ~  a 


(9) 


where  6f^  is  the  perturbed  distribution  function  of  the  a  species.  We 

assume  that  perturbed  quantities  vary  as  exp  [i  (k  •  x  -  u>t)]  where 

k=ke  +ke  +ke  and  k  =  -  i  3/3x  is  an  operator.  Equations  (4)  and 
^  xx  y  y  z  z  x 

(5)  becomes 


k2  £q>  =  4  ire  (6n^  -  Sn^) 


(10) 


k2 


6A  =  -  —  6J 
~  c.  ~e 


(11) 


In  writing  Eqs.  (10)  and  (11)  we  have  also  used  the  Coulomb  gauge 

(V  •  6A  =  0)  and  have  assumed  to2  «  c2k2.  Only  the  electrons  contribute  to 

perturbed  current  in  Eq.  (11)  since  u>2  «  c2k2. 

Pi 

We  now  calculate  6nQ  and  6Je  as  functions  of  the  perturbed  potentials 
using  linear  Vlasov  theory.  Details  of  this  calculation  are  given  in  the 
Appendix  A.  Making  use  of  5n^  and  6 Jg ,  we  write  Eqs.  (10)  and  (11)  as 


54>  +  Da  5A  +  Dx  6A  =0 
k^  4>x  x  <>z  z 


(12) 


k  k 

^  D  ,  4  +  D  5A  +  D  6A  *  0 
k  xij>  xx  x  k  xz  z 


(13) 


D  54  +f-  D  5A  +  D  5A  =  0 
z4>  k  zx  x  zz  z 

y 


(14) 
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In  order  to  solve  Eqs.  (12)  -  (14)  we  assume  that  k2  «  k2  and  expand 

x  y 

each  D  about  this  parameter.  That  is*  we  write 

D  =  +  (k  /k  )  +  (k2/k2)  Analysis  of  the  relative 

x  y  x  y 

magnitude  of  each  term  in  D  indicates  that  the  first  order  term  can  be 


neglected  when  u>/k  v  «  1.  Since  we  are  interested  in  modes  which 
b  ye 

have  (it  <  k  v.,  this  criteria  is  satisfied  when  T  /T,  »  (m  /a,)*^.  In  the 
y  i  e  i  e  i 

1/2 

opposite  limit  T  /T  «  (m  /m.)  ,  the  first  order  terms  cancel  exactly, 

e  i  e  l 

Thus,  we  can  write  D  =  +  (k2/k2)D^h  Making  use  of  this  relationship 

and  eliminating  and  SA^  from  Eq.  (12),  we  arrive  at  the  following 
second-order  differential  equation 


p(oj,ky  ,x) 


3  2  <5  <p 
3x2 


-  q((u,ky  ,x)k2  6<p  =  0 


(15) 


where  p  and  q  are  derived  and  defined  in  Appendix  B  (Eqs.  (B32)  and 

( B34) ) .  In  writing  Eq.  (15)  we  have  made  the  indentif ication  k2  +  3 2/3x2 

and  retained  terms  only  to  order  k2/k2.  Thus,  the  sixth-order  set  of 

x  y 

differential  equations  in  Eq.  (12)  -  (14)  is  reduced  to  a  second-order 

differential  equation.  The  crucial  assumption  in  this  analysis 


is  k2  «  k2. 
x  y 

Although  Eq.  (15)  is  very  complex,  in  general,  its  form  is  simple  and 
amenable  to  numerical  analysis.  Moreover,  in  certain  parameter  regimes, 
analytical  solutions  are  possible.  We  now  turn  our  attention  to  solution 
of  Eq.  (15).  We  first  present  an  analytical  analysis  which  highlights  the 
various  influences  of  finite  0  on  the  lower-hybrid-drift  instability  in  a 
sheared  magnetic  field.  We  then  present  numerical  results  for  a  broader 


parameter  regime 


lit.  Results 


A.  Analytical  Results 

Earlier  analytical  studies  ot  the  intluence  of  magnetic  shear  on 
the  lower-hybrid-drift  instability  in  a  finite  3  plasma  wave  were 
restricted  to  the  cold  electron  limit  (Tg  -*■  0  ) .  Recently  Huba  et  al. 
(1982)  1'*  extended  previous  analytical  results  to  include  7B  electron  drift- 
wave  resonances.  We  present  the  major  result  of  this  analysis  to  shed 
light  on  the  nature  of  the  two  important  finite  B  effects:  electromagnetic 
coupling  and  VB  drift-wave  resonances. 

Huba  et  al.  (1982)  derive  the  following  dispersion  equation  which 
describes  the  lower-hybrid-drift  instability  in  a  finite  B  plasma 
containing  a  sheared  magnetic  field  for  thi  lowest  order  mode 


k  V  w-k  V 

D(to,k)  =  1  +  k^p^  / 2 - +  i{/n  ( - + 

es  w  kVj 


kv,  up  T, 

i  pe  Kes  i  /  2\ 

- -  +  i  s  exp(-s  )A 

u  .  u  L  T  e  e  o 

pi  s  e 


(10) 


where 


J  (t  )  /Ui  (^  )1/Z  s  1/2 
W  kp  '  (2T~ }  se 
es  i 


Jl(cr> 


and  p2  =  p2  /(I  +  B./ 2),  p2  =  2T.  /m  SI2,  V,.  =  (v2/2S2J)  3  in  n/3x, 
es  Kes  i  es  iee  di  1  i 

2  2  1/2 

<d  *  Airne  /m  ,  Q  *  k  rT  s  ,  and  s  «  (w  /k  V..  )(2/B  )•  This  equation 
pa  a  r  y  Lc  e  e  r  y  di  y  e' 

2  2  2 

is  derived  based  upon  the  following  assumptions:  2u>pj/k  »  1, 

w  2  »  n  ?  s  »  1  w/k  v  »  1  k2r2  «  1,  V..  «  v.  ,  T  «  T. 
pe  e  e  II  e  y  Le  d  i  i  e  i 

and  Bj  «  1.  The  first  imaginary  term  in  Eq.  (16)  is  the  destabilizing 
term  due  to  Inverse  ion  Landau  damping.  The  second  imaginary  term  is  the 
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stabilizing  effect  of  magnetic  shear.  It's  origin  in  the  analysis  is  a 
term  «  (k  /co)2  in  the  magnetized  electron  response.  Physically,  magnetic 
shear  leads  to  stabilization  since  it  allows  wave  energy  to  propagate  away 
from  the  excitation  region  (i.e.  where  k ^  =  0).  The  final  term  is  a 
damping,  term  due  to  the  electron  VB  drift-wave  resonance 
The  real  frequency  is  given  by 


=  k  V  ,./(l  +  k2p2  / 2) 
y  di  es 


(17) 


where  shear  corrections  to  w  have  heen  neglected.  Tlie  mode  is  stabilized 

r 

when  Im  L)((u,k)  =  0  or 


T. 

/tt  t-  sfi  exp(-Se)A2). 
e 


where 

the 

subscript  refers 

to  the  critical 

value 

of  L^/Lg.  The  maximum 

value 

of 

the  RHS  of  eq. 

(18)  occurs  for 

k2p2 

es 

-  k2p2  =  2  so  that  all 
es 

wavenumbers  are  stable  when 


L 

f-"1 

1  L 1 


cr 


l+k2p2 


es 


/2 


(1  +  -±-) 

2k2p2 

es 


(ui 

V  V. 
1 


k2p2 

ejs 

l+k2p2  /2 
es 


(1H) 


(~) 

'•L  'cr 


4-  o  +  4"1)  tvd'1  +  lU  T1  se  e*P<-*fc>Aj] 
i  e 


(19) 


and  s  =  1/S  - 
e  e 

Two  interesting  points  concerning  Eqs.  (lb)  and  (19)  are  the 
following  First,  the  finite  S.  dependence  in  the  first  term  of  Eqs.  (18) 
and  (1S»)  arise  from  the  electromagnetic  correction  due  to  SA^  (i.e,  the 
transverse  magnetic  field  fluctuations)  The  influence  of  this  correction 
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is  to  increase  the  amount  ot  shear  necessary  to  stabilize  the  mode-  That 
is  as  0^  increases,  the  shear  length  Lg  necessary  for  stabi 1 ization 
decreases  so  that  the  mode  is  harder  to  stabilize  Physically,  this  occurs 
because  of  the  fluctuating  electric  field  associated  with  inhibits  free 

streaming  electron  flow  along  the  magnetic  field  which  reduces  the  rate  at 
which  energy  can  be  convected  away  from  the  localization  region.  Secondly 
the  final  term  is  Eq.  (19)  represents  the  resonant  7B  correction  which  is  a 
damping  effect.  This  term  tends  to  decrease  the  amount  of  shear  necessary 
to  stabilize  the  mode  Thus,  the  finite  6  corrections  have  different 

influences  on  the  shear  stabilization  criterion.  Lossely  speaking, 

electromagnetic  effects  are  destabilizing  (i  e.  a  stronger  shear  is  needed 
to  stabilize  the  mode  from  the  (5  =  0  situation)  while  the  resonant  7B 
effects  are  stabilizing  (i.e.,  a  weaker  shear  is  needed  to  stabilize  the 
mode  from  the  £5  =  0  situation) 
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B.  Numerical  Results 


We  now  present  numerical  solutions  of  Eq .  (15)  for  a  variety  of 

parameter  regimes.  We  solve  Eq  (15)  by  first  re-writing  it  as 

®  A*1  -q(a»,k,x)<S$  =  0  (2u) 

3x2  y 

where  0  =  q/p  and  use  a  finite  difference  scheme  (i.e  Numerov  method)^ 
to  obtain  eigenvalues  and  eigenfunctions.  The  boundary  conditions  used  are 

6$(x)  =  -\/a  '  t-xP  1"  i  J  dy  W^2(y)  J  for  lx|  >  “  (21) 

Q  '  (x) 

and  x6*  =  V  f°r  x  *  0. 

dX 

The  sign  of  the  WKb  solution  in  Eq .  (21)  is  chosen  such  that  a  damped 

solution  is  obtained  in  the  limit  [xl  ■*  <*>.  Hie  integrals  contained  in  Q 

are  performed  numerically  using  a  Chehychev  quadrature  method. 

Figure  2  is  a  plot  of  y/w..  vs  kp  for  V  /v.  =  1.0  T  /T.  =  1.0 
^  ih  es  di  i  e  i 

ui  /Ci  =  10.0  and  L  /L  =  0.1.  We  consider  two  values  ol 
pe  e  ns 

3  B  -  0.0  and  B  =  0.2  The  B  =  0. 0  curve  (dashed  curve)  is  the 

electrostatic  result  and  is  presented  as  a  reference  to  aid  in 

understanding  the  influence  of  finite  B  etfects.  Three  B  -  0.2  curves  are 

presented:  (1)6  *  0,  6ou  =  0;  (2)  6  =0,  *  0;  and 

r  em  VB  ’  em  Vb 

(3)  6  *  0,  6_,  *  0.  These  conditions  have  the  following  meaning.  In 

em  Vb 

calculating  Q  we  have  arbitrarily  included  two  coefficients  6  and  6„„ . 

em  VB 

The  parameter  6  modifies  the  electromagnetic  coupling  term  (5  w  /ck) 
em  em  pe 

and  the  parameter  5_„  modi! its  the  electron  Vb  drift  velocity  (6_„  V,  ). 

VB  VB  Be 

Bv  chosing  6  and  6„,.  equal  to  0  or  1,  we  can  understand  how  these 
'  °  em  VB  1 
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Plot  of  y/wou  vs.  kp  for  V. . /v.  **  1.0,  T  /T.  =  1.0,  u>  /J1  *  10.0 

Jlh  es  dt  1  ’  e  l  pe  e 

and  L  /L  ■  0.1.  Two  values  of  0  are  considered:  0  *  0.0  and  0  *  0.2. 
n  s 

The  parameters  6em  and  6^  refer  to  coefficients  that  modify  the 
electromagnetic  coupling  term  and  the  electron  VB  drift  velocity  (see 
text  for  detailed  explanation). 
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finite  0  effects  independently  influence  the  instability  in  the  presence  of 
magnetic  shear.  That  is,  by  setting  6  =  1  and  6^  =  0,  we  retain 

electromagnetic  coupling  effects,  but  neglect  electron  VB  orbit 
modification  effects.  Conversely,  by  setting  d  =0  and  dr,..  =  1,  we 
retain  electron  VB  orbit  modification  effects,  but  neglect  electromagnetic 
coupling  effects.  Neither  of  these  limits  are  self-consistent  but  are 
taken  for  pedagogical  purposes.  The  self-consistent  limit 

is  6  =1  and  &n..  =  1.  The  top  curve  considers  5  =1  and  d-,,,  =  0.  Note 

that  the  influence  of  electromagnetic  effects  is  to  increase  the  growth 
rate  relative  to  the  0  =  0.0  curve.  This  is  most  pronounced  in  the  long 
wavelength  regime  (kpgs  <  i)  since  w^/ck  is  largest  in  this  regime.  The 
bottom  curve  considers  d  =0  and  5,,,.  =  1.  Note  that  the  influence  ot  the 
electron  VB  drift  is  to  decrease  the  growth  rate  relative  to 


the  0  =  0.0 

curve . 

The  enhanced 

damping, 

due 

to 

the 

dissipative 

electron  VB 

drift-wave 

resonance,  is 

strongest 

i  n 

the 

short 

wavelength 

regime  (kpes>  1)*  This  is  because  the  perpendicular  resonant 

velocity  (v^f)  is  proportional  to  k  ‘  in  this  regime  and,  therefore,  more 

electrons  can  participate  in  the  resonance.  Finally,  the  self-consistent 

result  (5  =1  and  =  1)  lies  in  between  the  two  extreme  limits, 

em  VB 

Electron  VB  damping  dominates  in  the  short  wavelength  regime  while  there  is 
a  balance  of  the  electromagnetic  and  VB  effects  is  the  long  wavelength 
regime.  These  results  are  consistent  with  those  presented  in  Davidson  et 
al.  (1977). 2 

In  Fig.  3  we  plot  y/w„,  vs.  kp  for  V,./v.  =  l.U,  u  /Cl  =  1U.U, 

°  Jill  es  ui  i  pe  e 

L  /L  =  0.1,  3  =  0.2  and  two  values  of  T  /T- :  T  /T.  =  0.1  and  T /T,  = 

ns  eiei  ei 

1.0.  The  growth  rates  for  these  two  curves  are  comparable;  the  only 
significant  difference  occurs  in  the  short  wavelength  regime  where  the 
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l.U  than 


electron  VB  drift-wave  resonance  causes  weaker  growth  for  Te/Tj  = 

T„/T:  =  0.1.  This  is  due  to  the  tact  that  V„  <*  T  . 
e  i  Be  e 

In  Figs.  (4)  and  (5)  we  plot  the  wave  potential  Q  (curve  a)  and 

eigenfunctions  66  (curve  b)  for  V,./v.  =  1.0,  w  /Q  -  10.0,  L  /L  -  0.1, 

di  l  pe  e  ns 

3  =  u.2  and  kp^  =1.5  In  Fig.  (4)  we  consider  T^/Tj  =  0.1  and  in  Fig. 

(5)  we  take  T  /T^  =10.  The  eigenvalues  for  these  cases  are 

w/tdj^  =  0.83  +  1  0-16  for  T^/T.  =  0.1  (Fig.  4)  and  =0.58  +  1  0.15 

for  Te/Ti  =  1.0  (Fig.  5).  Although  the  growth  rates  for  these  two  cases 

are  comparable  (as  in  Fig  3),  examination  of  the  wave  potentials  ((<) 

indicates  an  important  difference  between  the  "small”  and  "large"  T  /Tj 

limits.  First,  in  Fig.  4a  we  note  that  Qr  possesses  an  "anti-well" 

character  for  x/p  <  4.0  while  0.  is  roughly  constant.  For  x/p  >  4.0 
es  I  es 

0r  begins  to  increase,  but  is  still  negative  when  6$  asymptotes  to  0 

at  x/p  ,  =  7.5  (Fig.  4b).  The  eventual  increase  in  Q  is  due  to  electron 

GS  * 

Landau  damping  which  would  not  occur  if  Tfc  -  0.  On  the  other  hand,  Qj 

increases  sharply  for  x/pes  >20  which  is  due  to  the  finite  growth 

rate-  Thus  Fig.  4  indicates  that  mode  localization  is  primarily  due  to 

the  outward  convection  of  energy  i rom  the  region  x  =  G  in  the  limit 

Tt>/T|  ♦  0.6  Stabilization  of  the  mode  can  occur  when  the  outward 

propagation  of  energy  is  faster  then  the  growtli  of  the  mode.  On  the  other 

hand,  for  hotter  electrons  (Fig.  5a),  the  "anti-well"  character  of  Qr  is 

only  weakly  evident  for  x/Pt,s  <  1.0;  for  X/Pes  >  1*0,  increases  and 

becomes  positive  for  x/p  >  2.2.  The  mode  is  localized  within  the 

es 

region  x/pgs  <  3.5  (Fig.  5b).  In  this  case  electron  Landau  damping  is 

playing  the  dominant  role  In  the  localization  of  the  mode.  As  the  mode 
propagates  outward  from  x  =  0,  it  is  rapidly  dissipated  locally  by  electron 
Landau  damping  Thus  Figs.  (4)  and  (5)  Indicate  that  different  processes 


2  4  6  8 

X/  £?es 

(b) 

Figure  4 

^L')t  os  tin;  wive  potential  aiul  eigenfunction  for  V  /v.  -  l.U, 

cl  I  i 

u  /Q  =  10.0,  L  /L  =  0.1,  0  =  0.2,  kp  =  1.5  and  T  /l.  =  U.l. 
pe  e  ’ns  es  e  i 

The  eigenvalue  Is  w/w^  =  0.83  +  i  0.16.  The  subscripts  r  and  i 

reler  to  real  and  imaginary,  respectively.  (a)  Wave  potential  Q.  (b) 
Eigenfunction  64). 


are  involved  in  mode  localization,  depending  on  the  value  of  Te/Ti4  either 

outward  energy  propagation  (Te  <<  Tj )  or  electron  Landau 

damping  (T^>  T\  )  localize  the  mode. 

In  Fig.  6  we  plot  Y/w^  vs.  Ln/Lg  for  V^/Vj  =  1.0,  ta  /Sl^  ~  10.0, 

kp  =  1  5,  3  =  d.O  and  1  0  and  T  /T.  =  0.1  and  1.0.  The  value  of  kp 
es  6  1  tiS 

cliosen  correponds  roughly  to  maximum  growth.  In  general,  the  curves 

indicate  that  y  <*  -  a(Ln/Lg)-rb,  wliere  the  slope  a  and  intercept  b  depend 

upon  3  and  T  /T . .  Two  important  features  of  this  curve  are  the 

following.  First,  for  Tg/Tj  =  1.0,  the  growth  rate  for  3  =  1  0  is  always 

less  then  that  for  3  =  00.  The  value  of  shear  necessary  to  stabilize  the 

mode  for  3  =  0.0  is  l_n /I-s  =  0.25  while  for  3  =  1.0  is  Ln/Lg  =  0.17. 

Second,  for  Te/Tj  =  0  1,  the  growth  rate  for  3  =  1.0  is  less  than  that 

for  3  =  0.0  when  I^/Lg  <  0.18,  but  is  greater  than  that  for  3  =  0.0  when 

L  /L  >  0.18  Thus  as  3  increases  in  the  low  electron  temperature  limit, 
n  s 

the  amount  of  shear  necessary  to  stabilize  the  mode  increases  (i.e.  Ln/Ls 
becomes  larger);  this  confirms  the  result  predicted  from  the  analytic 
theory  (Eq.  (19)). 

In  Fig.  7  we  plot  y/u)„,  vs .  3  for  V  /v.  =  1.0,  T  /T.  =  1.0, 

°  Ih  dii  ei 

a)  /ft  =  10.0,  kp  =  1.5  and  various  values  of  L  /L  .  The  main  result  is 
pe  e  es  ns 

that  as  En/Ls  increases  the  growtii  rate  y  decreases,  as  expected.  Also, 
as  Ln/Ls  increases,  the  slopes  of  the  curves  change,  with  a  plateau-like 
structure  developing  around  3  =  0.5. 

Finally,  we  comment  on  the  numerical  accuracy  of  the  results 
presented.  Solving  Eq.  (15)  requires  a  substantial  amount  of  computer 
time.  Typically,  100  -  200  grid  points  were  used  in  the  finite 

differencing  of  the  differential  equation.  At  each  grid  point  12  numerical 
integrations  are  required  to  obtain  Q.  And,  in  general  5-10  interations 
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ei-=? 
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*-s 

Figure  6 

Plot  of  y/«„,  vs.  L  /L  for  V,  ,/v.  =  1.0,  oj  />?  =  1(;.0,  kp  =  1.5, 

X.n  ns  dl  l  pe  e  es 

0  =  0.0  and  1.0,  and  T  /T.  =0.1  and  1.0. 

e  l 
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(3 

Figure  7 

Plot  of  y/«.u  vs.  3  for  V  /v  =  1.0,  T  /T  =  1.0,  to  /:7  =  ll» 

J6n  dl  1  e  i  pe  e 

kp  =  1.5  and  several  values  of  L_/L„. 
es  n  s 
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were  required  to  obtain  an  eigenvalue  tor  a  given  set  of  parameters.  Thus, 
approximately  10^  numerical  integrations  were  needed  to  obtain  a  single 
eigenvalue.  In  order  to  minimize  the  computer  time  used,  some  accuracy  was 
sacrificed  by  using  fewer  grid  points,  in  both  the  finite  difference  scheme 
and  the  numerical  integrations,  than  possible.  However,  higher-resolution 
grid  spacings  were  used  to  estimate  accuracy  for  different  parameters.  We 
find  that  tiie  qualitative  results  presented  are  reliable,  but  that  the 
quantitative  results  are  accurate  to  within  5  -  10%  of  the  actual  values. 
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IV.  Discussion 


A  self-consistent  theory  of  the  lower-hybr id-dri f t  instability  in  a 

finite  (3  plasma  containing  magnetic  shear  has  been  presented.  The  theory 

incorporates  the  important  finite  3  effects  of  (1)  elect rostat ic  and 

electromagnetic  coupling  and  (2)  electron  VB  drift-wave  resonances.  The 

main  conclusions  of  ttiis  study  are  as  follows.  First,  magnetic  shear  is  a 

strong  stabilizing  influence  on  the  lower-hybrid-drift  instability,  which 

is  a  well-know  result.^  Second,  finite  (3  effects  can  play  an  important 

role  in  determining  the  amount  of  shear  necessary  for  stabilization,  as 

shown  by  Eq.  (19).  Interestingly,  the  two  finite  3  effects  mentioned  above 

act  in  opposing  ways.  The  effect  of  electromagnetic  coupling  can  be  viewed 

as  destabilizing,  that  is,  it  tends  to  increase  the  amount  of  shear 

necessary  to  stabilize  the  mode  (see  Te/Tj  =  0.1  curves  in  Fig.  6). 

Physically,  this  occurs  because  electromagnetic  oscillations  generate  a 

fluctuating  electric  field  along  the  magnetic  field  (due  to  <SA||).  This 

inhibits  free  streaming  of  electrons  along  the  magnetic  field  which, 

therefore,  reduces  the  rate  at  which  energy  can  be  convected  away  from  the 

localization  region.  The  effect  of  the  electron  VB  drift-wave  resonance 

can  be  viewed  as  stabi  11  zing,  that  is  it  tends  to  decrease  the  amount  of 

shear  necessary  for  stabilization  of  the  mode.  This  is  illustrated  in  the 

Te/Tj  =  1.0  curves  of  Fig.  7  and  in  Fig.  8.  Physically,  this  occurs 

because  the  electron  VB  drift-wave  resonance  is  dissipative  and  reduces 

growth  In  the  localization  region.  The  key  parameter  which  dictates  which 

finite  8  effect  is  dominant  is  Te/Tj .  Electromagnetic  effects  dominate 

when  T,  <<  T< .  while  VB  resonance  effects  dominate  when  T  >  T. . 
e  l  e  l 

The  magnitude  of  Te/Ij  also  plays  a  role  in  the  nature  of  the  shear 
stabilization  of  the  lower-hybrid-drift  instability.  For  Te  <.<  Tj ,  the 
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Stabilization  occurs 


wave  potential  Qr  has  an  "anti-well"  character. ^ 

because  the  rate  outward  energy  propagation  (away  from  the  localization 

region)  exceeds  the  growth  rate.  On  the  other  hand,  for  T  >  T  the  wave 

e  1 

potential  Qr  has  a  "well"  structure  and  the  wave  is  dissipated  locally  by 
strong  Landau  damping.  This  is  shown  in  Figs.  (A)  and  (5). 

These  results  are  applicable  to  reversed  field  plasmas  whicli  contain 
magnetic  shear. ^  Perhaps  the  simplest  example  of  this  is  illustrated  by 
the  e<ju  i  J  i  br  i  um 


B  =  B  lee  +  tanh  (x/A)e  ] 
-  o  y  z 


(22) 


and 


n  =  n  sech2(x/A) 


(23) 


where  e  <<  1  Here,  the  z  component  of  the  magnetic  component  reverses 

direction  at  x  =  0,  but  the  total  field  remains  finite-  Studies  of  the 

lower-hybrid-drift  instability  in  reversed  field  plasmas  with  e  =  0 

indicate  the  following.  Based  upon  a  nonlocal  linear  theory,  “  the 

fundamental  mode  is  localized  away  from  the  neutral  line  (i  e.  ,  x  =  0)  at 

roughly  |x|  >  A  with  a  half -width  Ax  «  A.  higher  order  modes  also 

localize  about  |x|  =  A  but  have  a  much  broader  halt-width.  Ax  <  A. 

1/2 

However,  these  nodes  do  not  penetrate  chosen  than  |x|  -  X(T  /!'  )  ot  the 

p  el 

neutral  line  because  of  electron  VB  drift-wave  resonance  damping. 
Subsequent  work  on  the  evolution  of  a  reversed  field  plasma,  using  an 
anomalous  resistivity  model  based  upon  the  nonlocal  mode  structure  of  the 
lower-hybrid-drift  lnstabilfty,  found  that  magnetic  flux  is  transported 
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towards  the  neutral  line  and  that  the  current  increase  at  the  neutral 


20 

line.  Thus,  the  instability  can  eventually  penetrate  closer  to  the 

neutral  line  than  predicted  by  linear  theory. 

It  we  take  e  *  0,  then  based  on  the  definitions  of  1^  and  Ls  it  can  be 
shown  that 


L 

n 

L 

s 


_ e 

sinh( 2x/X) 


_ 1 _ _ 

e2  +  tanh2(x/X) 


(32) 


Based  on  Eq.  (31),  it  is  clear  that  IWL^  ♦  0  in  the  limit  x/X  ♦  «, 
while  L^/l^  «  when  x/X  >  0.  Moreover,  for  this  equilibrium 


3 - -1 - -1 -  (33) 

cosh2(x/X)  e2  +  tanh2(x/X) 

so  that  3  ♦  1/e2  >>  1  as  x/X  ♦  0.  Thus,  for  x/X  <  1,  both  magnetic  shear 
and  finite  3  are  important  to  the  stability  properties  of  the  lower-hybrid- 
drift  instability  in  a  sheared  reversed  field  plasma.  For  sufficiently 
strong  magnetic  shear  (i.e.,  e)  the  linear  mode  penetration 
distance,  |x|  ,  to  the  neutral  line  can  be  substantially  larger  than  when 
there  is  no  shear. ^  Thus,  magnetic  shear  can  inhibit  (or  perhaps  even 
prevent)  the  penetration  of  the  mode  to  the  neutral  line  even  in  the 
nonlinear  evolution  of  the  plasma.  This  can  be  significant  to  both 
laboratory  and  space  plasmas  where  anomalous  diffusion  in  the  field 
reversal  region  can  be  crucial  to  the  dynamics  of  the  plasma. 
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Appendix  A 


We  outline  the  derivation  of  the  perturbed  electron  density  (6n  )  and 

e 

the  perturbed  currents  (^Jxe  and  6Jze)  necessary  to  obtain  Eqs.  (12)  - 
(14).  We  first  calculate  the  perturbed  distribution  function  which  is 
given  by 


6fg  =  ~  /  dT  exP  [1  *  Jj  •  x'(t)  -  i(i)T ]  (-  i^k.  64>  +  6A 


l  3F 

+  ;(v'x(ikx  5A))  }  •  t~- 
C  "■  ">  3  v 


(Al) 


where  x'(t)  and  v'(t)  are  the  unperturbed  orbits.  That  is, 


^  A  A 

v ' ( t)  =  v.  cos(SJ  T-0)e  -  [v,  sin(8  t-0)  +  (v  ,2/2il  )3  An  B/3x)e  +v,e 

~  iexxe  1  e  y  l  z 


(A2) 

A  A 

x'(t)  *  (v^/Q^)  sin(neT-®)ez+l(vj_/^e)  cos(SleT-0)-T(v12/2»ie)3  An  Il/Sxje^ 

(A3) 

+  v  t  e 

I  z 

where  “  vx  +  Vy  and  “  e  B/mec.  The  unperturbed  electron  distribution 
function  is  chosen  to  be 

Feo  *  n  ^  ~  exp  [-v2/v2]  (A4) 

e 

where  v2  *  2T  /m  and  v2  ■  v2  +  v2  +  v2. 
e  e  e  x  y  z 

Substituting  Eqs.  (A2)  -  (A4)  into  (Al)  and  performing  the  temporal 
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integration  we  arrive  at 


6f 

e 


2e 

m  v2 
e  e 


[(«♦ 


ck 


K-  5Av  +  \T  6A J ) 


_  v 


)  exp  -  -r)(n-m)]  R  {J  J  6$ 

i.  e  n  m 

n,m 


(A5) 


u>-k  v  k  v  k2  , 

+  {  — j— -  Zr—J  J  +  —  —  (i  cos  \[>  J  J  -  —  J  J  sin  ill)}  6A 
ck  k  n  m  c  2  nmonm  r  x 

y  y 


v,  k  k 
1  x  z 


k  v  k2  , 

+  {(— j r  t~ - —  (l  +  — ))J  J  +  —  17-17“  U  cos  ♦  J  J  -  —  J  J  sini|/)}6A  ] 

ck  k  o'-  ,  2  n  m  ckk  rmn  onm  YJ  z1 

y  y  kp  y  y 


where  Rfi  =  (w-kyVde)/(u  -  n  fle~ky  VB(>(v2/v2)  -  kxvz),  t|<  =  tan_1(kx/ky ), 
the  argument  of  the  Bessel  functions  is 

°  =  klVl/fle’  Vde  =  "  (ve/2iie)3  *°  n/3x*  VBe  "  “  (^72^)3  *n  B/3x* 

and  v2  *  2T  /m  . 
e  e  e 

The  perturbed  electron  density  and  current  are  defined  as 


6n  *  / d3v  <5f 
e  J  e 


(A6) 


6J 


~e 


fd3v  v  6f 
*  ^  0 


(A7) 


Making  use  of  Eq.  (A5)  in  Eqs.  (A6)  and  (A7)  we  obtain 


29 


k.  k 

4V  -(2,e  »*  )  -*[(«♦  +  -J  +  -k‘“,)) 

y  y  y 


+A  /  d(u2)  exp  (-u2)  {J2  Z  6$ 
o 


o 

k  k  v  ,  ,  v  k2 

-  J2 
i  o  'ck 

y  y 


(A8) 


+  J2  (-£  Z  +  -  I  z'j-i  — -  --  COS  <|<  J  J.  c  Z)6A 

'■k  o  '•ck  kc2  1  ck2  o  1  e  ‘  x 

y  y 

+  (-£-  J2  Z+  Vi  (1^  )  i  z'-i  iiJijJ  co.  *CZ)«A  ] 
'•ck  o'-k  a>'-,o;2  ckk  ol  1  z 

y  y  ky  y  y 


00 

5J  =  -(2ne  )  *v  A  /  d(u2)  u  exp(*m2)  [-1  J  J  .cos  ty  Z  6<f>  (A9) 

xe  v  de J  e  4  o  1 


+  <-*  -  t*  2  *  ar  r  r4'Vi  c“  *  '  r  ^  J!  c°'2  K2} 

y  y  k 

77  y 


+  {_i  drV  z  +  "V’*1  +  3  I  z  )Vi  COB  ♦  “  c-  ir  IT  Ji  cos  ***><*«  1 
y  y  ky  y  y 

-i 

6J  =(2ne  X2  )  v  A  /  d(u2)  exp  (-u2)  [J2  Z  6<J> 

ze  de  e  J  o 

o 

k  k  v  v  k^  / 

+{-£-  r-5  J2  Z  <1  -  -*-i  c)-i  J  J.  —  —  cos  Il>  Z  }6A  (AID) 

ck  k  o  u  o  1  c  k2  x 

y 


y  y 


+  <dr  V  -  "V  +  ™  cos  ♦  Vi  ^rr  V>6azi 

y  y  ky  y  y 

where  A  *  (w-k  V,  )/k  v  ,  u  *■  v,/v  ,  the  argument  of  the  Bessel  functions 
y  de  z  e  i  e 

is  o  *  (k  ,  v  /!1  )u,  the  argument  of  the  Z  functions  is 
lee 

C  *  (u)-k  V„  u2)/k  v  and  X2  “  v2/2to2  •  In  writing  fcqs.  (A8)  -  (AID)  we 
*e  y  Be  z  e  de  e  pe 

have  retained  only  the  n  ■  0  terra  in  the  summation  over  cyclotron  harmonics 


since  u»  <<  Q  . 

e 
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The  perturbed  ion  density  in  an  unmagnetized  plasma  is  simply 


where  r. 

1 


n 

=  -  Y-  [1  +  Z(^±)l 

i 


=  (u>  -  kV  )/kv.  ,  V,  =  (v?/2Q.)3  An  n/3x  and  •  2T  /m  . 

ai  1  ai  l  l  ill 
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Appendix  B 


We  outline  the  derivation  of  Eq.  (15)  and  present  the  detailed  form  of 
p  and  q  used  in  the  analysis.  Substituting  Sn^  and  6Je,  derived  in  the 
Appendix  A,  into  the  Maxwell  equations  (Eqs.  (10)  and  (11))  we  obtain  the 
following  set  of  equations 


(Bl) 

(B2) 

(B3) 

(B4) 

(B5) 

(B6) 

(B7) 

(B8) 

(B9) 


and 


00 

K  -  A  /  d(u2)  exp(-u2)  J2  (k^r^u)  2  (Ce)  (*10) 

$9  0 

ao 

K.  =  A  /  d(u2)  exp(-u2)  u  JQ  (k^^u)  Jj  z  Ug  )  (bll) 

^  o 

00  i  ' 

K  =  A  /  d(u2)  exp(-u2)  J2  (k^^u)  Z  (Cg)  (B12) 

9  o 

oo 

K.  =  A  /  d(u2)  exp(-u2)  u2  J2  (k.t.  u)  Z  (c  )  (B13) 

xx  ‘  1  I  Le  e 

o 

00  l  ' 

K  =  A  /  d(u2)  exp(-u2)  u  J0(k_LrLeu^  Jl^kLrLeu^  2  2  CB 1 4 ) 

o 

“  , 

Kzz  =  A  /  d<u2)  <ixp(-u2)  J^k^^u)  "2  Cg  2  (ce)  •  (B13) 

o 


In  the  above, 


<e  = 


w-k  V„  u2 
y  Be 

k  v  ’ 
z  e 


w-k  V . 
y  de 

k  v 
z  e 


w-k  . 
y  di 

kv. 


;  v 


v?  }  tn  a 


di 


23, 


3x 


de 


2Q 


e  3  fcn  n .  v 

3x  5  Be 


e  3  In  b 


23 


3x 


-  V 
2  de 


u  ■  v./v  ,  A2  =  v2/3u2  ,  v2  *  2T  /m  ,  u>2  =  4w  ne  /m  ,  r,  ■  v  /SJ  ,  Jn  is 

1  e  de  e  pe  e  e  e  pe  e  Le  e  e  n 

/ 

the  Bessel  function  of  order  n,  Z  =  dZ/dc,  k2  =  k2  +  k2,  and  k  He  (x/L  ). 

l  x  y  z  y  s 

We  expand  D  in  the  small  parameter  k2/k2  and  find  that 


U  =  D(0)  +  (— X) 

-  -  -  V 

y 


(B16) 
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Dif  “  1 - -  —  T- 1  +  — - 

99  k2  X2  li  L  1  1  1  k2  X2  ^ 

y  de  y  de 

(B17) 

D<2)  =  D<2)  =  i  r - \ - K-^Z) 

<bz  z<b  k  X,  ck  d>x 

y  de  y  T 

(B18) 

D(2)  _  .  DU)  _  1  .  ^\e  k(2) 

4>z  z  4>  k  X,  ck  <t>x 

T  r  y  de  y  r 

(B19) 

D<»  .  1  -f-k  K<2> 

xx  2k2  XX 

y 

(B20) 

d<2>  -  -  d<2>  -  i  —t‘  k<2> 
xz  zx  2k2  xz 

y 

(B21) 

=  1  +  — pe 

ZZ  c2k2  22 

(B22) 

oo 

Kif  *  "  AkyrLe  /  d(u2)  exp(_u2)  u  JoJlZ(C  ) 

YY  }  0 

(B23) 

09 

Ki2)  *  A  1  kvrLe  ^  d(u2)  exp(-u2)u2[JoJrJ2)]  Z(Cg) 

v  7  o 

(B24) 

09 

*  -  A  k  r.  /  d(u2)  exp(-u2)uJ  J,  -i-  Z  (6  ) 

6z  y  Le  ;  r  o  1  2  e 

0 

(B25) 

00 

K^2)  -  A  k  ru  /  dU^expC-u^Jjj'zU^ 

(B26) 

"A  I  kyrLe  ^  d^u2 )exp(-u2 )u2 [ JqJ^J2 Z  (cg) 

(B27  > 

o 


34 


(B2B) 


v,(  2) 


"  AkyrLe  ^  d(u2)exp(-u2)uJoJ1  |  Ug) 


and  tile  argument  of  the  Bessel  functions  is  k  r^u. 


We  now  solve  Eqs.  (B2)  and  (B3)  for  6A  and  6A  in  terms  of  6<|>.  We 

X  2 


find  that 


<SA  =  -  -r~ 


k  D  D  +D 
y  $x  zz  xz  ijiz  . 
- 

D  D  +  D2 
XX  zz  xz 


5A  = 


D  D  -  D  Dx 
4>z  xx  xz  4>z 

D  D  +  D2 

XX  Z2  XX 


5<t> 


(B29) 


(B30) 


Substituting  (Bio)  into  (B29)  and  (B30),  and  then  substituting  (B29)  and 

(B30)  into  (Bl)  we  obtain,  to  lowest  order  in  k2/k2, 

x  y 


p(w,k  ,x)  — ' 5<t’ - q  (a),k  ,x)  k2  6<j>  =  0 

y  ax2  y  y 


(B31) 


where 


p  =  d  (r(0)/r*;0^  +  D^2)(r(0)/r^0)) 

4><t>  ij>x  v  nx  dx'  <f>z  '  nx  dx  ' 


(B32) 


-  d<“>  (r<2>  r‘°>  -  r(0>  r<2>)/r<u>  2 

4>x  k  nx  dx  nx  dx  J  dx 


+  d<°>  (r<2>  r‘°>  -  r<0)r'w  )/r<°>2 

tyZ  v  nz  dz  ri*  '  ri 


nz  dz  ;  dz 


,  -  d<°>  -  d<°>  (rc0)  /r<°>)  +  D<°>  (r'°>/r<°>) 

H  <p<p  <j>x  v  nx  dx  '  $z  '  nz  dz  ' 


(B33) 


and 
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(B34) 


AO)  B 

■  Di°> 

D(0) 

+  D<0> 

D<0) 

ax 

<frx 

ZZ 

xz 

4* 

.(2)  s 

-  b<2> 

D(0) 

+  l/2^ 

”i0) 

+  o<"> 

d<2)  ■ 

- 

ax 

$X 

ZZ 

xz 

l pz 

xz 

(j)Z 

XX 

(B35) 


r(0)  ,  r(0)  _  D(0)  D(0)  +  d(0)2 

dx  dz  xx  zz  xz 


(B36) 


r<2>  =  0(2)  D(°>  +  2  D(0>  D(2)  -  D(0>2  D(2)/D(0) 
dx  XX  zz  xz  xz  xz  zz  zz 


(B37) 


r<°)  =  D(0)  D(0)  -  D(0)  D(0) 
nz  $z  xx  xz  ox 


(B3b) 


r(2)  .  „(2)  0(0)  .  d(2)  d(0)_  d(0)  u(2)  +  D(0)  D(0)/D(0) 

nz  oz  xx  xz  <j>x  xz  <|>x  xz  ij>x  xx 


(B39) 


r(2)  =  (2)  D(0)  +  2D(0)  d(2)  _  u(0)2  d(2)  (0) 

dz  zz  xx  xz  xz  xz  xx  xx 


(B40) 


and  we  have  mode  the  identification  k2  =  -  32/3x2. 


